Magnetic field correlations in a random flow with strong steady shear 
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We analyze magnetic kinematic dynamo in a conducting fluid where the stationary shear flow 
is accompanied by relatively weak random velocity fluctuations. The diffusionless and diffusion 
regimes are described. The growth rates of the magnetic field moments are related to the statistical 
characteristics of the flow describing divergence of the Lagrangian trajectories. The magnetic field 
correlation functions are examined, we establish their growth rates and scaling behavior. General 
assertions are illustrated by explicit solution of the model where the velocity field is short-correlated 
in time. 

PACS numbers: 47.35.Tv, 47.65.-d, 94.05.Lk, 95.30.Qd 



I. INTRODUCTION 

The subject of our work is the magnetic dynamo that is 
the magnetic field generation by hydrodynamic motions 
in a conducting medium. We theoretically investigate the 
effect in a conducting fluid (plasma, electrolyte) where a 
random hydrodynamic flow is excited. The principal ex- 
ample of such flow is hydrodynamic turbulence (see, e.g., 
Refs. [HQ) responsible for the magnetic field genera- 
tion in different geophysical and astrophysical phenom- 
ena, see Refs. |3l-[l0|. We consider the case where the 
magnetic field grows from small initial fluctuations and 
examine the evolution stage when the magnetic field is 
weak enough, so that one can neglect feedback from the 
magnetic field to the flow. The stage where the flow 
is independent of the magnetic field is called kinematic. 
The kinematic approach becomes invalid when the in- 
creasing magnetic field begins to affect the fluid motions 
essentially. In this case the velocity field is strongly in- 
fluenced by the Lorentz force, so that the induction dy- 
namics is no longer linear. In most cases that leads to 
saturation of the magnetic field fluctuations maintained 
by the hydrodynamic flow. Though the magnetic field 
cannot be described by a linear equation in the regime, 
the kinematic stage produces magnetic structures simi- 
lar to those occurring at the saturation state, see, e.g., 
Rcf. [11]. A possible explanation of the fact is related 
to strong intermittency of the magnetic field that implies 
that the feedback is concentrated in restricted space re- 
gions where the magnetic field is anomalously strong. 

We assume that the random flow, exciting dynamo, 
is statistically homogeneous in space and time. Usually 
an additional assumption is made that the flow is statis- 
tically isotropic. If the velocity field is short-correlated 
in time then it is possible to derive closed equations for 
the magnetic induction correlation functions (see, e.g., 
Ref. [l2[). The pair correlation function for the case has 
been analyzed in Refs. [l3|,[l4[- The complete statistical 
description of the magnetic field for a short-correlated 
smooth statistically isotropic flow has been done in Ref. 
[l5j where growth rates and structure of spatial correla- 



tion functions were found. However, it is of interest to 
consider random flows with an average shear flow, that 
are widespread in astrophysical applications. Such flows 
are statistically anisotropic and need a special analysis. 
Here we examine the case where a steady shear flow is 
complemented by a relatively weak random component. 
We focus on the analysis of growth rates of moments of 
the magnetic field (magnetic induction), the degree of its 
anisotropy and on structure of the magnetic field correla- 
tion functions. We aim to relate the magnetic statistical 
characteristics to those of the flow, thus revealing the 
most universal features of the dynamo effect. The gen- 
eral assertions are illustrated by the model where the ve- 
locity field is short correlated in time, that can be solved 
analytically. 

An additional motivation for our research comes from 
dynamics of polymer solutions that in many respects is 
similar to magneto-hydrodynamics, see Refs. [l?], EH- 
Particularly, we have in mind the so-called coil-stretch 
transition [jl| (see also Refs. [13, IH|) that is an analog 
of the dynamo effect. A decade ago the elastic turbu- 
lence was discovered (see Refs. |23 - l24| ) that is a chaotic 
hydrodynamic motion of polymer solutions, the elastic 
turbulence can be realized even at small Reynolds num- 
bers as against the traditional hydrodynamic turbulence. 
The elastic turbulence is a natural frame for applying an 
extension of the dynamo theory to polymer solutions. 

The behavior of the magnetic field moments at the 
kinematic stage in the presence of the strong shear flow 
was established in Ref. [25} . However, to examine the 
spacial structure of the magnetic field, one should know 
its correlation functions, that are studied in the present 
work in the framework of the same general scheme as in 
Ref. 25j. To verify our general predictions we examine 
the analytically solvable model where a short-correlated 
in time random flow is excited on background of the 
strong stationary shear flow. 

The structure of our paper is as follows. In Section HTl 
we introduce basic relations needed to analyze the mag- 
netic field correlations and dynamics. We present general 
dynamic equation, give its formally exact solution and 
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discuss statistical properties of quantities entering this 
solution. In Section IIHI moments and correlation func- 
tions of the magnetic field are investigated. We relate 
its growth rates to the growth rates of the separation 
between two close fluid particles and establish principal 
space structure of the correlation functions. Section llVl is 
devoted to the model where fluctuating component of the 
flow is short correlated in time. We establish the growth 
rates for the model and analyze in detail the pair corre- 
lation function. The obtained results are in agreement 
with our general assertions. In Section [V] we outline our 
main results and discuss their possible applications and 
extensions. 



II. BASIC RELATIONS 

We consider the magnetic field evolution in a conduct- 
ing fluid (plasma or electrolyte) where hydrodynamic mo- 
tions are excited. Then the magnetic field dynamics is 
governed by the following equation (see, e.g., Ref. [26j ]) 

d t B = (B ■ V)v - (v ■ V)B + kV 2 B, (2.1) 

where B is the magnetic induction, v is the flow velocity 
and k is the magneto-diffusion coefficient, inversely pro- 
portional to the electrical conductivity of the medium. 
The flow is assumed to be incompressible, V • v = 0. 
We also assume that the magneto-diffusion term in Eq. 
(|2.ip is small in comparison with those related to the 
flow. We consider the case where the magnetic field is 
relatively weak and, therefore, its feedback to the flow is 
negligible. Then the relation (|2.1j) is a linear equation 
determining the magnetic field evolution in a prescribed 
velocity field, this regime is called kinematic. 

The hydrodynamic motion excited in the fluid is as- 
sumed to be random (turbulent) and the velocity statis- 
tics is assumed to be homogeneous in space and time. We 
examine the magnetic field growth from initial weak fluc- 
tuations distributed statistically homogeneously in space 
at the initial time t = 0. The correlation length of the 
initial fluctuations / is assumed to be smaller than the ve- 
locity correlation length 77. If we consider hydrodynamic 
turbulence then a role of the velocity correlation length is 
played by the Kolmogorov scale. At scales less than 77 the 
velocity field v can be treated as smooth. The magnetic 
growth (dynamo) can be characterized by moments of 
the magnetic induction that exponentially increase over 
time t: 

(|B(t)| 2n )ocexp(7 n t), (2.2) 

Here angular brackets mean averaging over space. The 
exponential laws (|2.2p are characteristic of the kinematic 
dynamo since the equation (|2.1[) is linear in the magnetic 
induction B in this case. 

One of our goal is to express the growth rates (incre- 
ments) 7„ in Eq. (|2.2p via statistical characteristics of the 
flow. The natural measure for the growth rates 7„ is the 



so-called Lyapunov exponent of the flow, A, that is equal 
to the average logarithmic divergence rate of close fluid 
particles. A special question concerns an n-dependence 
of 7„. If the magnetic induction statistics is Gaussian 
then 7„ cx n. Deviations from the linear law signal about 
intermittency of the magnetic field. The intermittency 
implies that high moments of the magnetic field are de- 
termined by rare strong fluctuations. 

There are two different regimes of the kinematic mag- 
netic field growth. The first regime is realized if all char- 
acteristic scales of the magnetic field are much larger than 
the magnetic diffusion length r ( i, r d — \J re/ A. The as- 
sumed smallness of the diffusion coefficient implies the 
inequality r\ ^> r d . We assume also I ^> r d , then the 
diffusion term in Eq. (12.11) is negligible at the first stage 
of the magnetic evolution that we call diffusionless. The 
magnetic force lines are deformed by the flow without 
reconnections in this regime. However, distortions of the 
magnetic field by the flow inevitably lead to producing 
scales of order r d in the field. After that the magnetic dif- 
fusion is switched on that admits reconnections. This sec- 
ond (diffusion) stage is characterized by the growth rates 
different from ones describing the diffusionless regime. 

Let us present a qualitative picture explaining the mag- 
netic field evolution at the kinematic stage. The initial 
magnetic field distribution in space can be thought as 
an ensemble of blobs of sizes ~ I. Then the blobs are 
distorted by the flow being stretched in one direction 
and compressed in another direction. In the isotropic 
case directions of the stretching and compression vary 
chaotically in space and time whereas in our case they 
are attached to the shear flow: the blobs are stretched 
mainly along the shear velocity and are compressed in 
the direction of the shear velocity gradient. At the first 
(diffusionless) stage the blobs are deformed without in- 
tersections and the magnetic field induction grows like a 
separation between close fluid particles since the equa- 
tion (|2.1I) at k — coincides with the equation for the 
separations. 

The diffusionless stage finishes when the characteris- 
tic blob width diminishes down to the diffusion length 
r c i- Then the diffusion is switched on that leads to two 
effects. First, the diffusion prevents further shrinking 
the blob widths, so that they remain of the order of r dl 
whereas the blobs continue to be stretched in the direc- 
tion of the shear velocity. Second, due to reconnections of 
the magnetic force lines admitted by diffusion the blobs 
start to overlap. As a result, new blobs of a characteristic 
longitudinal size ry are formed, see FigLU The magnetic 
induction in such blobs can be found by averaging the 
induction of a large number N of initial blobs, the num- 
ber ./V grows exponentially as time runs. Averaging over 
the large number of random variables leads to appearing 
an exponentially small factor ~ 1/yN in the amplitude 
of the magnetic induction. Besides, the amplitudes of 
the initial blobs remain to increase over time as a sepa- 
ration between fluid particles. We conclude that at this 
second (diffusive) stage the magnetic field is still growing 
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FIG. 1: Sketch of typical magnetic blobs during the diffusive 
kinematic stage. 



exponentially over time but slower than at the first stage. 

We consider the case where the steady shear con- 
stituent of the flow is much stronger than the random 
one. Quantitatively, the condition is written as the in- 
equality s A where s is the shear rate. Indeed, the 
Lyapunov exponent in a pure shear flow is zero, and its 
non-zero value is associated with presence of the rela- 
tively weak random constituent of the flow. The distorted 
magnetic blobs are elongated mainly along the shear ve- 
locity. However, they are tilted with respect to the ve- 
locity direction due to presence of the random velocity 
component, see Fig. Q] The tilt possesses the same dy- 
namics as the direction of the polymer stretching in the 
same flow, see Ref. [2fJ. Therefore the tilt angle <fi, see 
Fig. [U can be estimated as <j) ~ V s - The tilt angle de- 
termines typical ratio of the magnetic field components 
B y /B x ~ X/s <C f where the X-axis is directed along 
the shear velocity that varies along the F-axis. Thus, 
the ratio s/A characterizes an anisotropy degree of the 
magnetic field. 



A. Lagrangian dynamics 

To analyze moments and correlation functions of the 
magnetic induction we need a solution of the magneto- 
dynamic equation (12.11) for the induction field B(t) in 
terms of its initial value B, B = B(0). We use here a gen- 
eralization of the scheme proposed in Ref. [27| that ex- 
ploits the Lagrangian approach to the fluid motion. Pass- 
ing to the Lagrangian frame, that moves with the fluid 
particles, one excludes the advection term in Eq. (|2 . 1 1) . 
Next, instead of solving the equation with the Laplace 
operator V 2 it is convenient to introduce Langevin forces 
£ mimicing the magnetic diffusion. Then we pass to the 
following stochastic equation 



d t R = v(t,R)+£, 



(2.3) 



describing Lagrangian trajectories disturbed by the 
Langevin forces. The latter can be treated as white noise 



characterized by its pair correlation function 
(&(ti)&(*2)) =2KSij6(ti-t 2 ), 



(2.4) 



where k is the diffusion coefficient figuring in Eq. (|2.I[) . 

Rewriting the equation (|2.I[) in terms of the quantity 
R(t) one finds then the following formally exact solution 
ofEq. 



B{t,r)= W(t)B[R{0)]\ 



(2.5) 



where floors mean averaging over the ^-statistics, in ac- 
cordance with Eq. (I2.4[) . To find the quantity R(0) one 
should solve the equation (|2.3[) on the time interval (0, t) 
with the boundary condition R(t) = r, posed at the fi- 
nal time. By other words, one should track the magnetic 
field back in time along the disturbed Lagrangian trajec- 
tories. The matrix W(t) in Eq. (|2.5p is a chronologically 
ordered exponent 



W(t) = Tcxpl / dt't(t' 



(2.6) 



where E(i) is the velocity gradients matrix, Eji = diVj, 
taken at the time t and in the spacial point R(t). The ma- 
trix W, that we call an evolution matrix, can be treated 
as a solution of the equation d t W = EVF with the initial 
condition W(0) = 1. 

The evolution matrix W has some general properties 
that follow from the definition (|2.6I) . The determinant of 
the matrix W is equal to unity, since the velocity gra- 
dient matrix £ is traceless, tr£ = 0, that, in turn, is a 
consequence of the incompressibility condition V • v = 0. 
Therefore WiW 2 W 3 = I, where Wi,W 2 ,W 3 are eigen 
values of the matrix W. All the eigen values are positive 
since they are positive initially (W\ = W 2 = W3 = I 
at t = 0) and cannot turn to zero because of the same 
relation W\W 2 W 3 = 1. Let us order the eigen values 
in accordance with W\ > W 2 > W 3 , then W\ > 1 and 
W3 < 1. At times t > A -1 , we are interested in, typical 
values of In W\ , In W3 can be estimated as ±Xt, therefore 
W\ is exponentially large whereas W3 is exponentially 
small. An estimation for W 2 depends on details of the 
flow statistics. Any case, W\ ^> W 2 >• W3 at times 
t > A" 1 . 

In the framework of the proposed formalism, correla- 
tion functions of the magnetic field B have to be calcu- 
lated by averaging products of the factors (|2.5p taken at 
the respective points over the statistics of the noises £, 
besides averaging over space. Thus, say, the simultane- 
ous correlation function, 



F 9 



,r 2r . 



(B i (r 1 )...B j (r 2n )), (2.7) 



has to be calculated in two steps. First, one should sub- 
stitute the expression (|2.5[) into the right-hand side of 
Eq. (|2.7p and then average the resulting product over 
the ^-statistics given by Eq. (I2.4[) . this averaging catches 
diffusion effects. Let us underline that the fields £ have to 
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be treated as independent for all 2n factors in the prod- 
uct. Second, one should average the result over space. 
An averaging over scales < r\ (traced back to the initial 
time) gives statistics of the initial magnetic field fluctu- 
ations, and averaging over scales > r\ counts different 
realizations of S. Therefore the latter is equivalent to 
averaging over the velocity statistics. This logic was re- 
alized for the isotropic random flow in Ref. [la ]. 

In the diffusionless regime, realized at t <C 
A" 1 \n(l/rd), one can neglect diffusion effects. Then at 
calculating the moment (\B\ 2n ) one can take a product 
of the identical factors (|2.5|) where R is simply a La- 
grangian trajectory terminated at the point r at time 
t. Then \B(r)\ 2n w W? n \B\ 2n where B is taken at the 
origin of the Lagrangian trajectory. Here just the fac- 
tor W 2n is responsible for the exponential growth of the 
moments, and therefore, we can restrict ourselves to the 
estimation |J-B(r)| 2 "J ~ W 2n Bo n where Bo is the char- 
acteristic value of the initial magnetic field fluctuations. 
In the diffusion regime, realized at t ^> A -1 \n(l/rd), the 
situation is a bit more complicated. 

Let us first consider the second moment. Then we 
deal with two trajectories, R and _R', terminating at the 
same point r at time t, but characterized by independent 
noises £ and The second moment can be written as 
the following average 



((B) 2 ) = ( B T [R{Q)]W T W'B[R'{0)] 



(2- 



where the subscript T means transposition of the object. 
An appreciable contribution to the second moment is re- 
lated to the trajectories with \R(0) — R'(0)\ < I. Since 
|.R(0) - R'(0)\ < r) and \R(t) - R'(t)\ = 0, the difference 
AR = R— R' stays to be much less than r\ at any time 
from the interval (0,t) for such event. Then we obtain 
from Eq. (|2.3[) expanding the velocity up to linear in AR 
terms 



d t AR — EAi? + £ — £', 



(2.9) 



where 2 can be taken at the point R or R' ', no difference. 
A solution of Eq. (|2.9|) . equal to zero at time t' = t, is 
written as 

AR(t') = -W(t') J dtxW- x {ti)[i{t x )-i'{tx)]. (2.10) 

To calculate the second moment we should know the 
Ai?(0)-statistics. Since the separation AU(0) is a linear 
combination of £, £' it should be treated as a Gaussian 
variable at averaging over the ^-statistics, then its prob- 
ability distribution function is completely determined by 
the matrix M: 

Mij = LA J R l (0)A J R J (0)J = An J dh W^WT*, (2.11) 



the expression for M is derived from Eqs. (I2.4l2.10p . The 
matrix M is symmetric, and its eigen values are positive. 



Let us designate the eigen values as m\,m\,m\ and or- 
der the values as m\ > 777,2 > ma, the inequalities become 
strong, mi 3> mi 3> 7773, if At ^> 1. Let us stress that 
the directions of the eigen vectors of the matrix M are 
"frozen" at At ^ 1 |29l431| . Then the integral determin- 
ing mi is gained at t — t\ ~ A -1 and we arrive at the 
estimation mi ~ fdW 3 ~ ■ The integral determining 1713 
is gained at t\ ~ A^ 1 and therefore 7773 ~ r^. An estima- 
tion for m2 depends on the time dependence of Wi- If 
Wi increases then mi remains of the order of whereas 
it grows like m-2 ~ rjW^ 1 if W2 decreases. 

Now we find a probability that AH(0) is less than I 
in the diffusion regime, when t 3> A -1 \n(l/rd)- One can 
think in terms of the components of AJ?(0) in the basis 
attached to the eigen vectors of the matrix M . Since 
mi ^> I then a probability that the first component of 
AR(0) is less than / is estimated as I /mi ~ (l/TdjWz. If 
Wi grows as time runs then both, 771,2 and 7713, are of the 
order of Td and therefore a probability that the second 
and the third components of AJ?(0) are less than I, is 
close to unity. Then we find from Eq. (|2.8[) 



L|B| 2 J -BKl/rJWWz 1 , 



(2.12) 



where Bq is a characteristic value of the initial magnetic 
field fluctuations and we used the relation W1W2W3 = 1. 

The situation with the decreasing W2 is slightly differ- 
ent. In this case 7712 3> I at the diffusive stage and there 
appears an additional small probability that the second 
component AR(0) is less than I, the probability can be 
estimated as J/1712 ~ (l/i'd)W2- Then one would obtain 
(\B\ 2 ) - BlW^l/rd) 2 , instead of Eq. (l2~T2l . However, 
an integration over space (at the next step of averaging) 
kills the leading term due to the solenoidal nature of the 
magnetic field B. Therefore one has to take into account 
the next term in the probability distribution of Ai?2(0) 
that gives an extra small factor (l/m2) 2 in the probabil- 
ity. Thus, we arrive at 



[\B\ 2 \~B 2 (l/r d ) 4 W 1 W 2 



(2.13) 



Note that the expressions (|2.12I2.13[) are equivalent to 
ones obtained in the Fourier representation for the sta- 
tistically isotropic case in Ref. [151] . However, the ex- 
pressions ([2.1212. 13[) . written for real space, are correct 
for the anisotropic problem (we are investigating) as well, 
and are in fact more suitable for the problem. 

Let us turn to the high moments. One can recognize 
that a principal contribution to the average U-B| 2 ™J is 
produced by configurations where the 2n points R a (0) 
are divided into n pairs with separations < I in each 
pair. Because of the independence of the white noises 
£ a , the probability of such event can be estimated as a 
product of probabilities for the second moment, that is 



B\ 2n \ - l\B\ 2 \ n , 



(2.14) 



where the second moment is given by Eq. (|2.12|) or Eq. 
([2.13jl . We have ignored a combinatoric factor in Eq. 
([2.14jl being interested in the time dependence of the mo- 
ments. 
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B. Evolution matrix 

The next step in finding the magnetic field moments 
is averaging over the velocity statistics. Before doing it 
we should establish statistical properties of the evolution 
matrix (|2.6[) . Some universal properties of such matri- 
ces, that can be treated as products of a large number 
of random matrices, are well established, see Refs. [29i — 
HH, the properties are revealed at t ^> A -1 . However, we 
examine a strongly anisotropic case, with steady shear 
flow dominating. That forced us to modify the consider- 
ation of Ref. [15j where the isotropic case (statistically 
isotropic flow) was investigated. 

For the anisotropic problem, it is convenient to use the 
Gaussian decomposition of the evolution matrix W — 
T^ATr, where Tl and are triangle matrices, 





X 











i 


1 


X2 


, fa = \ Ci 


1 


i 







1 / 


U 


C3 





and A is a diagonal matrix. Since both triangle matrices, 
Tl and Tr, have unit determinants, the determinant of 
A is equal to unity as well. 

The matrices are written in the reference frame at- 
tached to the shear flow: the axis X is directed along the 
shear velocity and the axis Y is directed along the shear 
velocity gradient. Therefore the shear velocity is written 
as v x — sy where s is the shear rate. For our flow, which 
is composed of the steady shear flow and a random com- 
ponent, the matrix of the velocity gradients Sjj = is 
a sum of two terms related to the shear and the random 
components of the flow: 

T,ji(t) = s8j X Siy + (Tji(t). (2-16) 

The random matrix <jji is zero in average and should be 
characterized in terms of its correlation functions. The 
trace of the matrix is zero, tr<7 = (due to the flow 
incompressibility) . Let us remind that the Lyapunov ex- 
ponent A of a purely shear flow is equal to zero. Therefore 
A is sensitive to a though the random flow is weaker than 
the steady one. 

Substituting the decomposition W — T l ATr into the 
evolution equation dtW = T,W, one finds 

f^ l tf L = f^dtfL+dtAA-i+AdtTRT^A- 1 . (2.17) 

The terms in the right-hand side of Eq. (I2.17[) are the 
left-off-diagonal matrix, the diagonal matrix, and the 
right-off-diagonal matrix, accordingly. Therefore, one ob- 
tains a closed (non-linear) equation for the matrix Tl 
that leads to homogeneous in time statistics of the ma- 
trix. Next, we obtain for components of the diagonal ma- 
trix <9 t AA -1 expressions that are random variables with 
statistics homogeneous in time. Therefore InAi, lnA2, 
and In A3 (where A, are eigen values of the matrix A) 
are subjects of the central limit theorem. Typically, the 
variables behave linear in time t with coefficients of the 



order of A. The situation with the matrix Tr is slightly 
more complicated since there are the exponential factors 
in the last term of Eq. (|2.17[) . Therefore some compo- 
nents of the matrix Tr behave exponentially in time like 
the factors. 

Based on the leading role of the shear term in the ex- 
pression (|2.16[) . one obtains from Eq. (I2.17P for Tl a 
hierarchy \ 3> %\ S> Xi- Therefore in the principal ap- 
proximation in A/s the only component, cr yx , is relevant 
and the equation for the matrix Tl is reduced to a single 
equation for the component x- 

dtX = s - xV, (2.18) 

where a = a yx . We conclude that the variable x P os ~ 
sesses a homogeneous in time statistics, in accordance 
with our general expectations. Note that \ ~ s ll ^ 1 
as follows from Eq. (|2.18[) . Keeping the main in x contri- 
butions to the diagonal terms in Eq. (|2.17[) one obtains 
diag (pt AA -1 ) = (— x<r, X* 1 ! 0). Therefore in this approx- 
imation 

diagA^e-VM), dtP = XO- (2.19) 

If t > A -1 then typically p ~ At » 1. 

One concludes from the equations for d , £2 , C3 follow- 
ing from Eq. (|2.17[) that at At ^> 1 the variable £1 is 
"frozen" at a level of order unity whereas the variables 
C2jC3 grow exponentially, they can be estimated as e p . 
However, the combination (1O3 — C2, entering , is 
'frozen" at a level of order unity as well as £1. Based on 
the results obtained for the matrices Tl, A, and Tr one 
finds eigen values of the matrix W . In the main approx- 
imation in A/s we obtain W\ ~ e p , W% ~ 1, W3 ~ e~ p . 
The expressions together with Eq. (|2.19[) lead to the re- 
lation 

X=(d tP ) = (x<j), (2.20) 

where averaging is performed over the velocity statistics. 

Since the variable p is an integral over time of a ran- 
dom quantity with homogeneous in time statistics, see 
Eqs. (|2.18l2.19p . the p-statistics possesses some univer- 
sal features at At 3> 1. Namely, at the condition the 
probability distribution function (PDF) of p can be writ- 
ten in a self-similar form [321 ] 

P(p) cx exp[-tS(p/t)], (2.21) 

where S is the so-called Kramer function (or entropy 
function). The expression (|2.21[) is a manifestation of 
PDF's for the so-called intensive variables (see, e.g., [33| ) . 
The expression (|2.2ip implies that relative fluctuations of 
p diminish as t grows. 

Let us consider moments of the divergence of close La- 
grangian trajectories in our random flow. The equation 
governing a separation between the trajectories AR is 
dtARj — SjjAi?i, it can be obtained from Eq. (|2.9p by 
putting £ — > 0. A solution of the equation is AU(f) = 
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WAR(0). Therefore at t 3> A 1 we arrive at the estima- 
tion AR(t) ~ AR(0)e p . Then the moments of AR can 
be calculated in the saddle-point approximation (justified 
by the inequality Xt ^> f): 

(\AR\ n ) = I dp P(p)\AR\ n cx exp(A„i), (2.22) 

A„ = -S(ip„) + rwp n , where S'(tp n ) = n. (2.23) 

Thus, the exponents A„ are determined by statistical 
properties of the Lagrangian trajectories. Note that the 
Lyapunov exponent A can be formally expressed via A n 
as A = d\ n / dn\ n= Q. 

General statistical properties of the separation AR for 
the random flow with strong average shear were estab- 
lished in Ref. [2l|, in context of the single polymer 
dynamics in such flow. One expects a strong intermit- 
tency of AR(t), it should be revealed in faster than linear 
growth of A„ at large n, since the linear law A„ oc n is 
characteristic of Gaussian statistics of AR(t). 

III. CORRELATION FUNCTIONS 

To find a time dependence of the magnetic field mo- 
ments one has to perform an additional average of the 
expression (|2.14[) over space that is equivalent to averag- 
ing over the p statistics and statistics of initial magnetic 
fluctuations. Thus, the 2n-th moment of the magnetic 
field induction is written as 

(B 2n (t)) = J dp P(p) [B 2n (t)\ . (3.1) 

In our approximation, W\ ~ e p , W2 ~ 1 and, there- 
fore, in the diffusionless regime B(t) ~ e p Bo, whereas 
in the diffusion regime one obtains LI-^| 2 J ~ ^o(V r d) eP > 
as follows from Eq. (|2 . 1 2[) . Substituting the expressions 
into Eq. (|3.1[) and integrating over p (in the saddle-point 
approximation) we find j n = \2 n for the diffusionless 
regime and j n — X n for the diffusion regime. Thus, we 
related the dynamo increments (growth rates) introduced 
by Eq. (|2.2I) to the statistical properties of the flow. Our 
results can be summarized in terms of the estimations 

(\B(t)\ 2n ) ~ exp(X 2n t)B 2n , t < A" 1 ln(//r d ), (3.2) 
(\B(t)\ 2n ) ~ (l/r d ) n eM^nt)B 2 n , t > X~ x ln(l/rd). 

The main contribution to the moments (B 2n {i)) is as- 
sociated with the component B x of the magnetic induc- 
tion directed along the velocity of the shear flow, see 
Fig. [TJ Let us turn to moments of the component B y 
directed along the gradient of the shear flow, (B 2n ). The 
moments are much smaller than the moments ( x B 2n (t)), 
the smallness is caused by the strong shear flow. One 
finds from Eqs. (12.1512.191) that [Bl(t)\ = X 2 [B 2 (t)\. 
Thus, the variable \ is a measure of the magnetic field 
anisotropy, x _1 determines the tilt angle <fi of the mag- 
netic blobs to the shear velocity, see Fig. [T] Since the 



variable \ possesses homogeneous in time statistics, the 
factor x~ 2 does not produce a difference in the growth 
rates, that is the both moments, (B 2n ) and (B 2n ), are 
proportional to the same exponent exp("f n t). However, 
the prefactors at the exponents are different. To find the 
difference in the prefactors is not enough to know sta- 
tistical properties of the variable p that determine the 
exponent. Generally, one should know a mutual prob- 
ability distribution of the variables pit") and x(i) that 
is quite complicated object depending on details of the 
flow dynamics. However, one can establish an estimation 
for typical fluctuations x ~ S /A that follows from Eqs. 
(12.1812.201) . Therefore, say, |_-B*(*)J ~ (s 2 /A 2 ) (*)J ■ 

There is a question concerning moments of the third 
component of the magnetic induction, (B 2n ). To analyze 
their behavior one should take into account the compo- 
nents of the matrix Tl, that we ignored at investigating 
B x and B y . Then we conclude that the time dependence 
of {B 2n ) is characterized by the same exponents exp(7 n i) 
both at the diffusionless and diffusion stages. As to pref- 
actors, they depend on details of the flow statistics. 



A. Pair Correlations 

Let us consider the simultaneous magnetic field pair 
correlation function 

Fijfor) = {Bifan+^Bjfan)). (3.3) 

Here angular brackets mean, as previously, averaging 
over space (that is an integration over r% with the inverse 
volume as a factor). We assume statistical homogeneity 
in space of both the velocity and the initial magnetic field 
fluctuations, that is why the space average (|3.3p charac- 
terizes the magnetic field correlations in the whole vol- 
ume. We consider the case r -C 77, then we can use the 
smooth flow approximation. 

Again, we start from the presentation (|2.5|) . Then, 
analogously to the second moment, the pair corre- 
lation function (|3.3I) can be written as Fij (t) — 

(^[W lk B k [R(0))W^Bi[R'{0)}\y where the trajectories R 

and R' terminate at the points ri+r and t*i, respectively, 
at time t. Then one obtains 

Fij{t) = {\W tk {t)W 3 i{t)F kl [AR{{))]\) , (3.4) 

where Fij is the initial (at t = 0) pair correlation function 
of the magnetic field fluctuations and AR = R~R . The 
correlation length I of F is smaller than rj, that is why we 
can consider |AJ?| < rj. Then both evolution matrices in 
the expression Q3.4[) can be taken at the same point R. 
Averaging in Eq. (|3.4[) can be treated as averaging over 
the velocity statistics. 

The difference AR satisfies the same equation (|2.9[) 
provided |A_R| <C However, now we are interested 
in the solution with the final condition AR — r. Such 
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solution is written as 



properties. We define the Mellin transform as 



AR(t') = W{t')W- l {t)r 
~W(t') J dhW^ihMh) - {'(t!)], (3.5) 

instead of Eq. (|2.f Oj) . We immediately conclude from Eq. 
(|3.5p that the pair correlation function coincides with the 
second moment if r < r^. Therefore further we examine 
the case r ^> r c i where the second term in Eq. ()3.5|) is 
negligible and we find 



AR(0) = W'^ty. 



(3.6) 



To be more precise, the expression (13.61) is correct pro- 
vided y ~ r 3> rrf, that is implied below. 

There are two different regimes for the pair correlation 
function. If t < A -1 ln(Z/r) then AR(0) is typically less 
than I, the regime exists if r -C I- In this case the two La- 
grangian trajectories, R and R' remains typically within 
the correlation radius I at t — and the behavior of the 
expression (|3.4[) is insensitive to the separation r. There- 
fore the pair correlation function i*y practically coincides 
with the single-point average (BiBj) in this regime and 
, consequently, its time dependence is determined by the 
increment 7 = A2. 

If f > A -1 ln(//r) then AiJ(O) is typically larger than 
2 and only rare events where |AJ?(0)| < I contribute 
to the correlation function. Using the representation 
W = TlATr one obtains from Eq. (|3.6|l an expres- 
sion AR(0) e p (r x — x r y)j where and r y are co- 
ordinates of the separation r. A probability that the 
quantity is less or of the order of I is estimated as e~ p l/r 
(provided r x ~ r y ~ r), that is an interval of values of 
X where Ai?(0) < I. Therefore BiBj ~ B^efl/r and, 
consequently, ~ Bp exp(Ait)Z/r. 



Let us collect the obtained results 



F(t)~Bgexp(A 2 f), 
F(t) ~ Bl exp(Ait)//r, 



f < A-Mn(//r) 
t > A _1 ln(Z/r), 



(3.7) 



where the inequality <C r -C 2 is assumed. Thus, the 
pair correlation function is governed by the same expo- 
nents as the second moment. In addition, we find an 
r-dependence of the pair correlation function. Note that 
the expression (|3.7p turns into the expression (|3.2[) for 
the second moment at r ~ as it should be. 

Returning to the expression p.4[) we conclude that a 
difference between the pair correlation function F^ and 
the moments (BiBj) is solely in the behavior of AR. 
Therefore relations between components of Fij controlled 
by the evolution matrices in Eq. (I3.4[) is the same as for 
the moments (BiBj), say, F yy ~ (A/s) 2 F xx . 



F(t,k) 



dr 
r 



— ik 



F(t,r), 



(3.8) 



where a direction of the radius-vector r is implied to be 
fixed. Due to smoothness of the velocity field, differ- 
ent harmonics F(t, k) evolve independently, being rep- 
resented as a sum of exponents characterizing different 
structures of F^. At times t ^> A -1 only the principal 
exponent survives, that is F(k,tp) cx exp[7(fc)f]. In order 
to return back to the real space one should perform the 
inverse Mellin transform 



F{l. r) = j T^exp 



— ik In - 

r 



l{k)t 



P(k). (3.9) 



The quantity F(k) is determined by the initial magnetic 
field fluctuations, correlated on the scale I. That is why 
we incorporated the quantity into the relations (|3.8I3.9[) . 

Some words about analytical properties of fi(k). We 
assume that at r > I the initial pair correlation function 
J-(r) decreases fast as r grows. Then the integral (|3.8[) 
converges if Im k > 0. Therefore F(k) is analytical in 
the upper fc-semiplane. Besides, the integral diverges (at 
small r) as k — > 0. Therefore singularities of T(k) lie 
in the lower fc-semiplane, starting from the point fc = 
0. A character of the singularities depends on analytical 
properties of the initial function T(r). If it is analytical 
in r then one expects that .F(fc) has a series of poles 
along the lower imaginary semiaxis, and the first one lies 
at fc = 0. Note that in accordance with general rules the 
integration contour in Eq. (I3.9[) should go above the first 
singular point fc = 0. 

Now we can draw some general conclusions taking into 
account that 7(fc) ~ A. If ln(Z/r) > Xt then the integral 
(|3.9[) is determined by a narrow vicinity of the point fc — 
0. Then F(t, r) cx exp[7(0)£] and we identify 7(0) and A2. 
If ln(Z/r) < Xt then the integral (|3.9[) can be calculated in 
the saddle-point approximation. To find the saddle point 
one should shift the integration contour into the upper 
semiplane to reach the saddle-point fc = iq+, where q* 
determines the minimal value of j(iq) at q > 0. Indeed, 
the increment j(iq) is real and therefore the point fc = 
iq+ is a solution of the equation d'y/dk = giving an 
extremum of the exponent in Eq. p.9[) . Then F(t,r) cx 
(l/r) q * exp[7(i(7*)i], and we identify j(iq*) with Ai, see 
Eq. (|3.7p . Note that in accordance with the asymptotic 
law (I3.7P should be equal to unity, q* = 1. 



B. Mellin transform 

It is instructive to examine the Mellin transform of the 
pair correlation function, the analysis reveals its scaling 



C. High Order Correlations 

Here we consider higher order correlation functions 
of the magnetic field (|2.7p . One obtains expressions 
like Eq. (|3.5p for separations Ai? between the points 
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Ri (0) , . . . , R 2n (0) that are needed to calculate F 2n in 
accordance with Eq. (l2"3j) . If At < ln(Z/|Ar|) for all 
separations between the points ri , . . . , r 2n then all sepa- 
rations AH(0) are less than In this situation we arrive 
at the same expression F 2n ~ B 2n (e 2np ) as for the 2n-th 
moment, and we conclude that F 2n oc exp(A2 n t), see Eq. 



Now we turn to the case where At > ln(l/r) where 
all separations Ar are assumed to be of the same order. 
Let us consider first the geometry where all the points 
7*1, ... , r 2n lie on a line, that is all vectors r a — rp have 
the same directions, and one can write Ar ~ r, where 
r is one of the separations. Then one arrives at the es- 
timation |Ai?| ~ e p {r x — X r y)i similar to the pair corre- 
lation function. Thus, one obtains the same probability 
~ e~ p l/r that the separations |AH(0)| < I. Then we 
find F 2n ~ Bl n {e i - 2n -^ p )l/r, where the factor e 2np origi- 
nates from the product of the matrices W, appearing in 
accordance with the expression (|2.5p . Averaging the ex- 
pression, one obtains F 2n [t) <~ Bq 71 exp(X 2n -it)l / r . Let 
us stress that the growth rates here are different from the 
growth rates of the corresponding moments. 

However, the above expression is correct only if t < 
A -1 ln(i/fd). For larger t the diffusion contributions to 
the differences A a R become relevant, see Eq. (|3.5I) . 
Then, manipulating by Xi it is possible to make less 
than I only one difference among AiJ(O). After that 
all other differences acquire typically values of the or- 
der of rd,e p , and a probability that a difference is smaller 
than I is estimated as (l/rd)e~ p . Thus we conclude 
F 2n ~ ((iMil/r^e^) cx exp(A„t). 

The same results (upto combinatoric factors) are true 
for a collinear geometry where the set r\, . . . , r 2n is sep- 
arated into n pairs with parallel vectors r a — rp char- 
acterizing the pairs. Then the corresponding differences 
AR behave like previously and the same arguments are 
correct. Let summarize our results for the collinear ge- 
ometry 

F 2 »(t) ~ B 2n exp(A 2n t), t < A" 1 ln(//r); 

F 2n (t) ~ B 2 n exp(X 2n -it)l/r, (3.10) 
A _1 ln(f/r) < * < A" 1 ln(Z/r d ); 
F 2n (t) ~ B^il/r^l/uY- 1 exp(A„t), * > A" 1 \n{l/r d ). 

It is a generalization of the expressions f|3 . T[) for the pair 
correlation function. Note also, that, as for the pair 
correlation function, in the collinear geometry the esti- 
mate By ~ (X/s)B x determines relations between dif- 
ferent components of F 2n whereas z-components need a 
separate investigation. 

If the collinear geometry is destroyed, then it is impos- 
sible to put all the separations |AU(0)| inside the scale 
I at any \, if t > A -1 ln(Z/r). In this situation a behav- 
ior of the correlation function F 2n is non-universal being 
sensitive to details of the initial space distribution of B. 
Any case, a value of F 2n in the non-collinear geometry 
is much less than in the collinear one. The situation re- 
sembles one realized for the randomly advected passive 



scalar on scales larger than the pumping length [36|. We 
conclude that at t > A -1 \n(l/r) correlations of the mag- 
netic field are concentrated near collinear geometries, de- 
caying away from the geometries. The decaying length 
is estimated as le~ xt at t < X~ 1 \n(l/r c [) and as at 
t> X-HnQ/rd). 



IV. SHORT CORRELATED FLOW 

Here we consider a strong steady shear flow comple- 
mented by a random component shortly correlated in 
time. The case admits an analytical solution and can be, 
therefore, used to varify our general assertions and pre- 
dictions. In addition, the case is naturally realized since 
the strong shear destroys correlations of the random com- 
ponent, that is why we expect that the short-correlated 
case is widely spread in real flows. 

In the short-correlated case the matrix of the velocity 
gradients a describing the random component of the flow 
has to be treated as white noise, that is a variable 5- 
correlated in time. In the isotropic case one obtains the 
following tensorial structure 

{v%k(ti)<7jn(Pd)=D(45ij5i m —5ik5j n —6i n 5jk)5(ti—t2), (4.1) 

where the factor D characterizes strength of the random 
flow and the numerical factor is introduced as in Ref. 
35]. However, as we already argued, the only relevant 
component of the random velocity gradients matrix in 
the case A <C s (that is a manifestation of the random 
flow weakness) is a = a yx . We characterize its statistical 
properties by the expression 



(a{ti)a(t 2 )) = 4DS(t!-t 2 ), 



(4.2) 



formally coinciding with Eq. (|4.1[) for the yx-component. 
Other components of a can have correlation functions dif- 
ferent from the expression (|4.1[) . The random component 
can be treated as weaker than the steady shear flow pro- 
vided D <C s. 

Statistical properties of the separation between close 
Lagrangian trajectories AR(t) in the random smooth 
flow with strong shear component in the short-correlated 
case are investigated in Ref. [37] (in the context of poly- 
mer dynamics). We present here results obtained in the 
paper without a derivation. Note that our variable x 
is related to the tilt angle (f> of Ref. 137] as x — c °t 4>i 
X ~ in the case of small tilt angles, we are interested 
in. An expression for the Lyapunov exponent found in 
Ref. [33 is 



A = 



r(i/6) 



£>l/3 s 2/3_ 



(4.3) 



Thus, the condition s ^> D guarantees the inequality 
A <C s, indeed. Note also that A ^> D and that A — » 
as D — > 0. The last property is a natural consequence of 
zero value of the Lyapunov exponent for a purely shear 
flow. 
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For the short-correlated case, it is possible to find the 
exponents A„ characterizing growth rates of the moments 
of AR(t), see Eq. ([2722]) . if n > 1. Then a saddle-point 
(instanton) approximation in the functional space can be 
used, see Ref. [34j, that leads to 



7n 



2 5/3 



n 4/ 3£) l/3 a 2/3 „ An 4/3 



(4.4) 



see Ref. [37j. The non-linear dependence of A„ on n, 
A„ oc n 4 / 3 , signals about strong intermittency of the flow. 
Note, that in our anisotropic case the growth rates A„ in- 
crease as n grows slower then in the isotropic case, where 
for the short-correlated flow A n oc n 2 , see Ref. jl5j |. 



A. Pair Correlation Function 

Next, we examine the two-point simultaneous corre- 
lation function (I3.3[) for the short-correlated in time ve- 
locity field. In the case it is possible to derive a closed 
equation for the correlation function (see, e.g., Ref. fl2|). 
We study here the function at scales that are much larger 
than the diffusion scale (but much smaller than the ve- 
locity correlation length rf) . Then it is possible to neglect 
diffusion effects and we omit all terms with the noise £ 
in subsequent relations. 

Let us shortly explain a derivation of the equation. 
First of all, it follows from the definition (12.61) that 



W(t) = Texp(H)W(i-r), E 



d£ % (4.5) 



where r is an arbitrary time (less than t). We choose 
t to be much smaller than A -1 but much larger than 
the velocity correlation time (it is possible for the short- 
correlated flow). Then H is a small factor, though the 
two factors in Eq. (|4.5p can be treated as statistically 
independent. Substituting the expression (|4.5p into Eq. 
(|3.4p , expanding the result into a series over S (upto the 
second order) and averaging the result inside the time 
interval (t — r,t) in accordance with Eq. (I4.2[) . one ob- 
tains a variation of F^j at passing from t — r to t. Since 
the variation is small, it can be rewritten in terms of a 
differential equation. 

Assuming the isotropic correlation function of the fluc- 
tuations (|4.ip one obtains the following equation 



dtFij = -2syd x Fij 



+4D 



SijFkk — —Fi 



rkdjF lk - r k diF] 



-(rV)i^ + ir 2 V 2 F„- - -r m r n d m d n F i:j 



(4.6) 



In the absence of the shear (at s = 0) the system of 
equations (14.61) leads to a closed equation for the trace of 
the correlation function H = F k k'- 



The equation coincides with one presented in Ref. [14[ 
(for the case of the scaling exponent £ = 2 and zero 
forcing) . 

Now we rule out irrelevant terms in Eq. (|4.6p using the 
following properties: s is much larger than D, character- 
istic value of x is much larger then that of y and, respec- 
tively, d y ^>d x . Then we can keep in Eq. (|4.6p solely the 
terms originating from a yx . The resulting equation leads 
to a closed system of equations for the three components, 
F XX i Fxy, F yv , of the pair correlation function: 

d t F xx = -2syd x F xx + 2sF xy + ADx 2 d 2 F xx , (4.7) 
dfF X y — '2syO x F X y ~\~ sFyy -\- c ZF)x dyF X y AF)xOyF xx , 

d t Fyy = -2 S yd X Fyy + 2DX 2 d 2 y Fyy - 8DxdyF X y + DF XX . 

Further, we use the dimensionless time T — (8Ds 2 ) 1 ^ 3 t 
and introduce the designations / = F xx , g — 
(s/Dy/ 3 F xy and h = {s/D) 2 / 3 F yy . 

Let us investigate a special case of coinciding points. 
At r = all terms with derivatives drop from the equa- 
tions (|4.7p and they take the form 



1 
1/2 
2 




(4.8) 



A growing solution of the solution is 




(4.9) 



The behavior corresponds to the increment 7 = 
(815s 2 ) 1 / 3 of the magnetic field second moment. From 
the other hand, the case corresponds to small r that is 
to the condition hx(l/r) > Xt. Therefore 7 = A2, and we 
conclude that A2 = (815s 2 ) 1 / 3 in our case. 

It is convenient to pass into the "polar coordinates" 
g, tp in the shear plane: x = gcostp, y = {D / s) 1 ^ gsvcup. 
Then we perform the Mellin transform f,g,h —> f,g,h 
in terms of g and derive the equations for f,g,h from 
the system (|4.7p . In terms of the quantity q — —ik the 
equations are written as 

&rf = f" cos 4 <p + f [sin 2 tp — 2(q + 1) sin ip cos 3 <p\ 
+f [q sin ip cos ip — q cos 4 p + (q 2 + q) sin 2 p cos 2 p\ + g, 
d t g = g" cos 4 tp + g' [sin 2 tp — 2(q + 1) sin tp cos 3 tp] + 
+g [q sin tp cos tp — q cos 4 p + (q 2 + q) sin 2 p cos 2 tp\ 

o -. ~ 1- 

— 2 cos tp f + 2q sin tp costp /+ —h, 

dji = h" cos 4 <p + h' [sin 2 tp — 2(q + 1) sin ip cos 3 tp\ + 
+h [q sin p cos p — q cos 4 p + (q 2 + q) sin 2 p cos 2 tp] 

-4 cos 2 tp g' + 4g since cos g + 2f, (4.10) 



d t H = D (10H + Qrd r H + r 2 d 2 H 



where prime designates a derivative over the angle tp. 
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FIG. 2: Increment of the Mellin transform of the pair corre- 
lation function on the imaginary axis. 



Next, we study the time evolution of the system (|4.10p 
numerically for different real values of q (so that k = iq is 
imaginary) using implicit difference scheme on the inter- 
val (— 7r/2, 7r/2) for ip with the periodic boundary condi- 
tions. We have chosen as initial conditions for /, <?, h the 
same Gaussian functions centered near cp = and with 
width of order unity. Then we extract the principal in- 
crement "/(iq) — (8I?s 2 ) 1 / 3 c((?) dominating the behavior 
of the system at T > 1. The dimensionless quantity c 
was extracted as 

1 /(Tp + T) 

c= T ln -m^T' 

where Tq + T is chosen to be large enough (near 30) 
and To is introduced to exclude an influence of an initial 
transient process, we have chosen To = T. 

The quantity c is plotted as a function of q in Fig. [2l It 
turned out to be positive everywhere, having a minimum 
at q = 1, c(l) = 0.435. The value q = 1 is in accordance 
with Eq. (|3.7|) and the general analysis given in Section 
IIIII As we argued there, the minimum value of y(iq) 
determines Ax, that is Ax = c(q^)(8I?s 2 ) 1 / 3 . The value 
of c at q = is c = 1, in accordance with Eq. (14.9[) 
and general arguments given in Section IIIII Thus, the 
obtained results confirm our general assertions. 



V. DISCUSSION 

We have analyzed the kinematic dynamo stage when 
the small-scale fluctuations of the magnetic field grow in 
a steady shear flow complemented by relatively weak ran- 
dom velocity fluctuations. The weakness is characterized 
by the inequality s ^> A where s is the shear rate and A is 
the Lyapunov exponent of the flow. Universal features we 
are established are revealed at times t A -1 . The shear 
makes the flow strongly anisotropic, which, paradoxically, 
simplifies the analysis of the dynamo phenomenon since 



a single component of the random velocity gradient ap- 
pears to be relevant. We analyzed the situation where 
the correlation length / of initial magnetic field fluctua- 
tions is less than the velocity correlation length r\ (that is 
the Kolmogorov length for developed turbulence). Prob- 
ably, the smallness of I is not crucial for our scheme since 
small scales of the magnetic field distribution in space 
are inevitably produced by the hydrodynamic motion. 

Let us underline that in the main approximation in 
A/s our problem is reduced to a purely two-dimensional 
velocity field (with components along the shear velocity 
and along its gradient). We have proved an existence of 
the dynamo in this case (that is the exponential growth 
of the magnetic field moments). The result obviously 
contradicts the statement of Refs. [38M40l ] (Zeldovich 
theorem) that there cannot be the magnetic dynamo in 
two-dimensional flows. We assert that this statement is 
wrong and the error is in ignoring the third component 
£?3 of the magnetic induction (perpendicular to the ve- 
locity plane). The third component satisfies the passive 
scalar equation and, consequently, decays exponentially. 
However, B3 cannot be ignored in the divergent less con- 
dition V.B = since the characteristic scale of the mag- 
netic field along direction of its growth increases faster 
than the magnetic field itself. One can check that all the 
terms in VB = d x B x + d y B y + d z B z decay with the same 
exponent, and therefore the condition "VB = leads to 
an effectively divergent in-plane magnetic field. The dy- 
namo effect is not forbidden for such field. A detailed 
analysis of the discrepancy will be published elsewhere. 
An existence of the dynamo effect for two-dimensional 
flows is a subject of numerical verification. 

For our general conclusions, we do not specify statis- 
tical properties of the random flow exploring only its 
smoothness at scales less than the velocity correlation 
length 77. Then it is possible to relate the kinematic 
growth rates 7„ of the magnetic field, see Eq. (|2.2I) . to 
intrinsic characteristics of the flow characterizing diver- 
gence of close Lagrangian trajectories, see Eq. (|2.22l) . 
We find that j n — \2 n in the diffusionless regime and 
7n = A„ in the diffusion regime. We also related the 
anisotropy degree of the magnetic field to the same intrin- 
sic characteristics of the flow. Thus, the main features of 
the magnetic field statistics (including its intermittency) 
are dictated by the flow statistics. Note that our general 
scheme can be applied without essential modifications to 
the statistically isotropic flows or to random flows with 
other types of anisotropy. 

We established main features of the magnetic field 
correlation functions. The pair correlation function be- 
haves like the second moment at small separations r, 
and increase with the growth rate, characteristic of the 
diffusion regime, at larger r, it is proportional to 1/r 
there. As to higher order correlation functions, the sit- 
uation is more complicated. As small time t they be- 
have as corresponding moments. However, at larger time, 
t > A -1 ln(^/r), correlations are peaked near the collinear 
geometry (where 2n points are separated into n pairs 
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with parallel separations) and there is an intermediate 
asymptotic regime when the correlation functions grow 
with rates that do not coincide with the growth rates of 
the moments. Then, at times t > A -1 ln(Z/r,j) the cor- 
relation function grows with the same exponent, as the 
corresponding moment in the diffusion regime. The scal- 
ing behavior of the correlation functions in the collinear 
regime is oc 1/r. Correlations decay fast at deviation 
from the collinear geometry. That reflects a complicated 
spacial structure of the magnetic field that is strongly 
correlated for special geometries produced by affine ge- 
ometrical transformations from initial magnetic fluctua- 
tions. 

Our general assertions can be verified by solving the 
model with short correlated in time fluctuating compo- 
nent, that admits a series of analytical results. The non- 
linear n-dependence of the growth rates 7„, j n oc n 4 ' 3 , 
at large n signals about strong intcrmittency of the mag- 
netic field. Therefore only rare events concentrated in a 
restricted fraction of space contribute to high moments 



of the magnetic field. We analyzed in detail the pair 
correlation function of the magnetic field for the model, 
the analysis confirms all our general assertions, including 
scaling behavior in different regimes. 

The ideology and the analytical approach developed in 
our work can be expanded to the dynamics of polymer 
solutions possessing elasticity that is described similarly 
to the magnetic field. Moving along these lines, we hope 
to clarify some aspects of the so-called elastic turbulence 
[2214241 ] that are still not explained. 
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